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Abstract 

A master equation description of the inherent dynamics is used to calculate the frequency- 
dependent linear thermo-visco-elastic response functions of a glass-forming liquid. Prom the imag- 
inary parts of the isobaric specific heat, isothermal bulk modulus, and isobaric thermal expansion 
coefficient, we define a quantity At p (uj) with the property that At p (u>) = 1 is equivalent to having a 
one-parameter description of the linear thermo-visco-elastic response. This provides an alternative 
to the well-known criterion based on the Prigogine-Defay ratio. 
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[Preliminary version] 

Several articles have examined if it is possible to give a one parameter description of 
the thermodynamic relaxation processes of viscous liquids The so 

called Prigogine-Defay ratio tests for a one parameter description by comparing 

the responses of the liquid state with the responses of a glassy state. This approach is 
problematic because the responses of the glassy state depends on the history of the glass. 
In other words the original Prigogine-Defay ratio is not an state function. Later tests 
compare the equilibrium (zero frequency] responses of a liquid to the instantaneous (infinite 
frequency) responses of the same liquid p,|f]]. In principle this should lead to an test quantity 
which is a state function. However instantaneous responses very difficult to measure, as one 
has to deal with experimental problems which are related to heat flow and thermal stresses 
in the sample. 

In the following we propose a new test which allows testing for a one parameter hypothesis 
by measuring the imaginary parts of three complex response functions at a single frequency. 
This way we can completely avoid the measuring of an instantaneous response. We examine 
a harmonic linear perturbation of a thermoelastic system with temperature To, pressure po, 
entropy Sq, and volume V . Writing p = p + Re(6p e lu)t ) and T = T + Re(5T e luJt ) the 
oscillation of V and S is described by V = V + Re(SV e iult ) and S = S + Re(SS e iwt ). 
The complex coefficients 5p, 8V,8S, and 5T are related by a matrix of complex, frequency 
dependent response functions. 
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V a p (uj) V k t (u) 

Here c p (u) denotes the isobaric specific heat capacity per volume, a p (u) denotes isobaric 
thermal expansion coefficient, whereas kt{oj) denotes the isothermal bulk compressibility. 
The matrix symmetric due to an Onsager relation, which will be proven in the course of this 
article. 



MARKOVIAN INHERENT DYNAMICS 



We shall base the following calculation one the so called energy landscape picture of the 
dynamics of a viscous liquid 8] . We assume that at a sufficiently low temperature a small 
molecular system spends most of its time vibrating around a minimum of the potential 
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energy function. Furthermore, we assume that the rate of transitions is so low that a system 
achieves vibrational equilibrium between two transitions. An ensemble of such systems is 
well described by a vector of probabilities P = (Pi, P2..., Pat) where P n denotes the fraction 
of the ensemble which is vibrating around the energy minimum n. For a large system the 
rate of transitions is proportional to the system size, and therefore we cannot expect a large 
system to achieve vibrational equilibrium between each transition. However assuming that 
a transition from one potential minimum to another minimum is a localized event, leads to 
an assertion of local vibrational equilibrium. 

In the following we examine a linear experiment in which we control the intensive variables 
T and p and calculate the resulting change of the extensive variables S and V. Here S denotes 
the ensemble entropy whereas V denotes the ensemble average of the volume. 

S = -E P ™(^ + ^ lnP -) ( 2 ) 

m 

G n (T,p) denotes the vibrational Gibbs energy of a single minimum, whereas the Gibbs 
energy of an ensemble is given by 

G(T,p,P) = J2 P n(G n (T,p) + k B TlnP n ). (4) 

n 

We introduce the notation X = (T, —p) and Q = (S, V), and note that the following equation 
is fulfilled for all choices of P 

dG 

Q a = -q^, a = 1,2. (5) 
The equilibrium distribution is easily found by minimizing G under the constraint Yl m — 

Wp) = exp (_^) . ( 6 ) 

For a linearly perturbed system, the rate of entropy change, S, is related to the heat flow, 
J q , by the relation J q = TS j^j. Thus the rate of change of the energy is given by 
U = TS — pV. We insert in the relation G = U — TS + pV to get 

G = -Sf + Vp = -^Q Q X Q . (7) 



When the thermal energy per degree of freedom is low compared the height of an energy 
barrier, the transitions between minima is well described as a Markov process, and the 
dynamics of P m it is described by a master equation P|. 

p = X^w p 

rn 

The rate matrix W depends on T and p, and the equilibrium distribution P^} fulfills 
J2m W nm Pm = 0. We proceed to examine linear perturbations around a reference state 
X° = (T°, -p°). Furthermore we introduce W^ n = W mn {T ,p ) and P^ = P%>{T ,p ). We 
let 5 denote perturbations around the reference state to get 

5P n = Y, W nJP m + 5W nm Pl 

m 

We note that 5(J2 m W nrn P%) = 0. This leads to £ m <W nm P£ = -EmC^- We 
insert to get 

5P n = W nm (SP m - 5P%) . (8) 



A little algebra and the use of © leads to [2J 



d In P eq 

m 



(dQl 



Qf3+5 P1 k B y (9) 



dX p k B T \dP, 
We insert in (jHJ) and remember Ylim ^n m Pm — to get 

SPn = £ Wl m \5P m - £ ^T P ^ 6X ^j ■ ( 10 ) 
The steady state fulfills 

« - iuSln) SPn = £ j^WLP^SXp. (11) 

n m,f3 

Here 5P n and 5Xp denote complex amplitudes. For uj > we can simplify (jll)) by defining 

Anmico) = -4- £ (W° ~ iul)-J W? m Pl (12) 

Here / denotes the identity matrix. We substitute (fT^j) into (|TT|) to get 

SP n = J^AnU^^rSXp- (13) 
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This equation describes the dynamics of the ensemble. In order to calculate the response 
functions, we define = — dx d2G dXfj and expand the complex amplitudes 6Q a to the first 
order to get 

= E i£r 5p » + E ( 14 ) 

n n (3 

We insert in ()13j) to get an expression for the thermoviscoelastic response functions that 
were introduced in ((TJ). This gives SQ a = J a p{uS)8Xp, where 

JafiM = J% + E l§? A ™^7§r- ^ 

m,n 

As A nTO (oo) = we get J a p(oc) = J^%. We can show that J a p(uj) is symmetric by introducing 
a symmetric matrix Y mn with the elements Y mn = (P^)'^W mn (P^)^ Q. This gives 



Anm(co) = -^-^(P^M + y) n ^ m (P:)i (16) 

The matrices (iul + Y) 1 and Y commute so it is easy to show that A nm (u>) is symmetric. 
We insert in (|15|) to conclude that the matrix J(u) is also symmetric. 

PRIGOGINE-DEFAY RATIO 



Havingcalculated the response functions we can insert in the "linear Prigogine-Defay 



ratio'' 



n = (cp(0) ~ Cp(oo))(kt(0) - « r (oo)) ^ 



T (a p (0) - a p (oo)) 2 
We insert (JT5Jl and remember A nm (oo) = to get 



= 7 as av^ ' ' 



From the Schwartz inequality to find II > 1. We can describe the special case II = 1 by 
noting that null space of A mn {ui) is given by null(A(o;)) = span{e}. Here e is a vector 
with elements on the form e n = ■h 2M- We n °te that II = 1 is only fulfilled we can define 



constants 7 and c, so that all the following equation is fulfilled for all n 



9V dS 

9R=^R +C ' (19) 



This result is similar to earlier treatments of the Prigogine-Defay ratio. As mentioned in 
the introduction, the Prigogine-Defay ratio is very difficult to measure, and therefore we 
propose the following ratio which is inspired by 



very diff 



A Tp (u) 



KM) 2 ■ 



We insert (JT3J) to get 



A Tp (u) 



ds A ii ( \ as 

>m,n dP n ^nmK^J Qp m 



/V 9V A" (:.:) dV 
yZ^,m,n dP n nm\ r 1 J dP m 



dS 



A': 



m,n dP n nm 



[CO. 



8V 
dP m 



From the Schwartz inequality we find that Atp{&) > 1, for all u g]0, oo[. In order to analyze 
the reverse relation we note that null(A"(ci;)) = span{e}. This tells us for all u> g]0, oo[ that 
A Tp (u) = 1 is equivalent to (fIT?|) and thus A Tp {oj) = 1 is equivalent to II = 1. We can get a 
better understanding of the implications of if we approximate G n by a Taylor expansion 
in T to get G n = H n — TS™ b . This approach is essentially equivalent to a series of recent 



articles on thermodynamics of the energy landscape 



141 ]. A little calculation yields 



dS H n — G — k B T 



(20) 



We insert in (fT^j) to get 



dV 



Jo 
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G + k B T 
T 



(21) 



This formula shows that if IT = 1 and if G n = H n — TS™ , then the equilibrium fluctuations 
of inherent enthalpy are directly correlated to the fluctuations of inherent volume. 



SINGLE PARAMETER MODEL 



In the following we will describe a simple formalism that makes it easy to perform cal- 



culations on a model fulfilling II = 1. 
formalism for thermodynamics 3> ^6 



we introduce the following parameter 

5 s 



Jl7 Jlc 1 

mm 



his calculation is inspired by the circuit model 



18(. Inspired by Goldstein and La Nave et. al. 
x dV 
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OR 



-5P m - 



(22) 



We combine (|22|) with (|14|) to find that in a linear experiment the instantaneous values of 
S and V are functions of T,p and fe. 



SS = de + J^5T - J£5p 
SV = 7 fe + J£5T - J™5p 



(23) 



This model is different from the order parameter models introduced by Meixner and Pri 



101 ] . In the following we will 



gogine, as it allows for 5e to have non exponential relaxation 
show that the dynamics of this system can be described as the flow of Gibbs energy between 
the two reservoirs 1 and 2 with the Gibbs energy densities G\ and Gi- These reservoirs will 
be chosen so that the total Gibbs energy G is given by 



G — G\ + G*2 



(24) 



We define the function G\ to be the loss Gibbs energy during an isostructural process leading 
the system, from a state (T,p, P) to the state (T Q ,p , P). 



Gt 



T ,p,P pT,p,P 

V dp* - S dT* 

T ,Po,P JT ,p,P 



(25) 



We insert ()23|) and = to find G\ is a function of ST, Sp and 5e. Furthermore we get 
the following second order expansion. 

G!(5T,5p,5e) = V 5p - S 5T 



I 


5T ^ 


T 


' TOO 
J ll 


TOO 

J 12 


l\ 


/ 


5T 




—dp 




TOO 

J 21 


TOO 

J 22 
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V 


5e j 
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V 


5e 



J 



We would like to derive an equation of conservation of Gibbs energy for reservoir 1. Therefore 
we define V = and note from (J23} that S = and V = fff to get 



Gt = —ST + Vp — ipe. 
A similar equation for reservoir 2 follows from (|7j).()24|) and (|26j) 

G 2 = G-G 1 = i/je 



(26) 



(27) 



Most of the literature on network thermodynamics is expressed in terms of energy flows, and 
in order to stay consistent with this literature we define reservoir energies U\ = Gi + ST—pV , 



and U2 = G2 to get 



u = U l + U 2 

ir 1 = TS-pV-i)e 

U 2 = 



(28) 



These formulas describe the dynamics of the liquid as energy flow between two reservoirs. 
We proceed to find equations for the dynamics of each reservoir. The dynamics of the first 
reservoir follows from (|25|) . We use the definitions of S, V and i/j and examine a harmonic 
perturbation to get 



TOO TOO 1 



ss 

SV 
-8$ 



TOO TOO 

J 2l d 22 1 



The partial antisymmetry of this matrix reflects the fact that —p, T and — ip are generalized 
forces on reservoir 2, whereas V, S and e are generalized charges. In order to describe the 
second reservoir, we need a relation between Sip and Ss. We use ([lip. (|12|) and (fT5)l to get 
Ss = a(u)5ijj, where (2(| 

dS 

or, ^dK, 



ST 
—Sp 

Ss 



(29) 



a(uj) 



(30) 



Finally we insert Se = a{uj)Sip in (J29j) to find a simple expressions for the three complex 
response functions that we defined in (jT|). 

Vo , v 



J™ + 7 2 a( W ) 
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